Macroscopic approximation to relativistic kinetic theory from a nonlinear closure 
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We use a macroscopic description of a system of relativistic particles based on adding a nonequilib- 
rium tensor to the usual hydrodynamic variables. The nonequilibrium tensor is linked to relativistic 
kinetic theory through a nonlinear closure suggested by the Entropy Production Principle; the 
evolution equation is obtained by the method of moments, and together with energy-momentum 
conservation closes the system. Transport coefficients are chosen to reproduce second order fluid 
dynamics if gradients are small. We compare the resulting formalism to exact solutions of Boltz- 
mann's equation in 0+1 dimensions and show that it tracks kinetic theory better than second order 
fluid dynamics. 
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I. INTRODUCTION 

Fluid dynamics successfully describes the long- 
wavelength, low-frequency dynamics of various systems 
including the quark-gluon plasma (QGP) created in ul- 
trarelativistic heavy ion collisions [lH8(. It is usually 
linked to the more fundamental kinetic theory descrip- 
tion by giving a closure, i.e., an expression for the one 
particle distribution function / in terms of the hydro- 
dynamic variables. The most widely used closures are 
those provided by the Chapman-Enskog expansion and 
the method of moments due to Grad (see e.g. 0), later 
on adap ted to relativistic systems by Israel and Stewart 
[if! These two approaches rely on writing / as a 

function of the same variables that describe the system 
at equilibrium, namely energy density e and velocity u 11 
for a system of relativistic massless particles as consid- 
ered here. 

Traditional fluid dynamics derived from kinetic theory 
by these methods has two important limitations. First, 
it relies on an expansion in gradients of hydrodynamic 
variables, which necessarily implies that the system is 
sufficiently close to equilibrium so that these gradients 
are small. In turn, this means that the system is also very 
close to being isotropic in momentum space. Second, it 
breaks down at large shear viscosity-to-entropy ratio rj/s. 
Comparison to relativistic kinetic transport calculations 
[H, [H| shows that the applicability of Israel-Stewart (IS) 
theory to describe heavy ion collisions becomes marginal 
at rj/s > 0.2 — 0.3, especially at very early times. 

The application of fluid dynamics to relativistic heavy 
ion collisions shows that the assumption of small gradi- 
ents breaks down at the earliest times after impact, and 
that large momentum-space anisotropies actually survive 
for many fm/c. Besides, near the edges of the fireball cre- 
ated in these collisions the momentum-space anisotropies 
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are large at all times (T3] ■ These conclusions are also sup- 
ported by studies based on the Anti de Sitter/Conformal 
Field Theory (AdS/CFT) correspondence, which show 
that viscous hydrodynamical behavior is reached at times 
when the system still possesses large momentum-space 
anisotropies, and that these anisotropies persist during 
the entire evolution of the plasma |l5l [HI . 

With respect to the second limitation of fluid dynam- 
ics mentioned above, one should bear in mind that what 
is actually extracted from data by comparison with vis- 
cous fluid dynamic simulations that use a temperature- 
independent rj/s is an average value of this ratio. In real- 
ity, T)/s depends on the local temperature of the medium, 
dropping from high temperatures to a minimum at the 
critical temperature, and risin g w ith decreasing tempera- 
ture in the hadronic phase [l7l1l8j. In this phase, rj/s can 
reach values ~ 0.5 — 1.0, which are too large for standard 
fluid dynamics to be reliable. Beyond this value of rj/s, 
an appropriate treatment of the dynamics of this matter 
should switch smoothly from a viscous hydrodynamic de- 
scription to a kinetic one 0-0, EH H3] ■ 

The purpose of this paper is to present a nonequilib- 
rium effective theory (NET) capable of handling highly 
nonequilibrium situations and large values of t]/s. The 
model was developed by the authors in [2l[ and applied 
to heavy ion collisions i n 12211 . The link with kinetic the- 
ory was established in [231 ]. and the inclusion of color 
degrees of freedom was carried out in [24[ (we note that 
in this paper we do not consider colored particles). 

Here, we want to focus on the comparison of the NET 
to second order fluid dynamics of |25H27| (we will call 
this theory SOFD, and note that it includes IS theory 
as a special case) and to kinetic theory. From numerical 
simulations of 0+1 boost invariant flow, we show that the 
NET tracks kinetic theory better than SOFD, even in the 
presence of very large momentum-space anisotropies or 
large values of 77/ s. Due to these features, the model can 
be useful to study in a unified way the dynamics of matter 
created in heavy ion collisions, from the early stage where 
large momentum-space anisotropies occur, through the 
deconfinement transition, beyond which large values of 
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r\l s arise if the temperature dependence of r\ and s is 
taken into account, until freeze out. 

In the context of heavy ion collisions, there have been 
previous extensions of standard fluid dynamics aiming 
at alleviating the above mentioned limitations (in other 
contexts there have been numerous extensions of stan- 
dard fluid dynamics, see Section [TTJ) . In |2gj . an extension 
of IS formalism to third order in hydrodynamic gradi- 
ents was carried out, and an approximate resummation 
to all orders was presented as well. It was found that the 
resummed theory tracked kinetic theory better that the 
third order formalism and much better than IS, even for 
(very) large values of r\j s ~ 3. In [29j], the IS equations 
were rederived in a novel way based on the exact def- 
inition of the dissipative part of the energy-momentum 
tensor instead of on the truncated expression used by IS 
[lol Hlj . This resulted in the IS evolution equations but 
with different expressions for the transport coefficients 
that greatly improved the agreement of this formalism 
with kinetic theory, even for large values of rj/s ~ 3 
as well. In |30l - t32| . the method of moments was ap- 
plied to an anisotropic distribution function instead of 
an isotropic one as it is usually done when deriving hy- 
drodynamics from kinetic theory. The result was a new 
formalism termed anisotropic hydrodynamics (AH) which 
can handle highly nonequilibrium situations and large 
values of rj/s. We note that our NET is closer in spirit 
to AH than to the developments of [28| and [2!| , because 
both the NET and AH include a nonhydrodynamic vari- 
able in addition to the usual hydrodynamic ones. As it is 
discussed later, this nonhydrodynamic variable accounts 
for the backreaction of the distribution function on the 
evolution of the hydrodynamic variables. 

This paper is organized as follows. In Section [TT] we 
describe the entropy principle on which the NET is 
based, review the closure obtained from this principle 
and present the evolution equations. In Section IIIII we 
compare the results of numerical simulations of the NET, 
SOFD and kinetic theory. This section contains our main 
results. In Section [IVI we end up with a brief summary. 



II. NONEQUILIBRIUM EFFECTIVE THEORY 

In this section we briefly review the NET that was 
developed in dlHHlj (see also [13). 

The NET is obtained from the closure picked up by 
the so-called Entropy Production Principle (EPP) (a re- 
cent review can be found in [33[ ) . This closure is not tied 
up to the usual hydrodynamic gradient expansion, and 
it is nonlinear in a nonhydrodynamic variable ^ v . The 
latter models the backreaction of / (that could describe 
a highly nonequilibrium situation) on the hydrodynamic 
modes (that relax much more slowly) . In the EPP, 7 MI/ is 
identified with the Lagrange multiplier of the variational 
problem whose solution gives the / that extremizes the 
entropy production. Physically, one can think of the EPP 
as selecting the dynamics of 7 A " / (accounting for the back- 



reaction of / on (e, u^)) in such a way as to extremize the 
production of entropy during the evolution of the system. 
This interpretation for 7 is rooted in nonequilibrium sta- 
tistical physics, particularly in the study of fluctuations 
around stationary states [3Jt36|. 

The idea of a closure based on a variational princi- 
ple whose Lagrange multipliers are macroscopic variables 
finds application in the formulation of Extended Ther- 
modynamics [37| . which use as a variational principle 
the Maximum Entropy Principle (MEP) instead of the 
EPP. On the other hand, the introduction of nonhydro- 
dynamic variables into a macroscopic description derived 
from kinetic theory is not unique to the EPP. Related 
approaches are discussed in [38|-|4Cj for the MEP, in 41 1 
where the moments of the collision term of the Boltz- 
mann equation are treated as independent variables (see 
also [i^]), in (30l - [32l | where AH is derived from an ex- 
pansion of / around an anisotropic state, in |43| in the 
framework of divergence- type theories (see also |2l| - [23l |). 
in [44| where / is expanded in moments to all orders 
in Knundsen number, and in [45j in the context of fluid 
models derived from the AdS/CFT correspondence. 

We note that various effective theories derived from 
the MEP and the EPP have found applications in the 
study of electron transport through semiconductors 13^- 
lifll . |4(| and radiative heat transport in a photon gas [42j . 
It was found that these effective theories are able to re- 
produce kinetic theory results better than standard hy- 
drodynamic models. 



A. From kinetic theory to the NET 

We consider a relativistic kinetic theory of massless 
particles in flat space-time with signature (—,+,+,+), 
and for simplicity neglect quantum statistics. The Boltz- 
mann equations is 



p^f(x,p).^=I col [f] 



(1) 



where / is the one-particle distribution function, / co i[/] 
is the collision kernel, and a semicolon stands for a co- 
variant derivative. We write the distribution function as 



/ = /o[l + (l + /o)xWo[l+x] 



(2) 



where /o = exp[— u u p u /T]; the four velocity is obtained 
from the energy momentum tensor through the Landau- 
Lifshitz prescription u^T^" — e(T)u v (T is the tempera- 
ture). We also define 



(...) 



Dpf (---) 



(3) 



where Dp = d 4 p8(p 2 )/(2ir) 3 . We shall base our develop- 
ments on a relaxation time approximation for 7 co i, suit- 
ably generalized to fulfill momentum conservation. We 
write the collision operator as 



1 

2t, 



R[FR X ] 



(4) 
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where r r is the relaxation time, F = F(u>) is an arbitrary 
function of energy ui = —p^u^, and R is an operator 
enforcing energy-momentum conservation (p^/coi) = 0. 
This form for / co i guarantees that the Second Law holds 
exactly. Moreover, it includes the well-known kinetic 
models of Marie (which is a relativistic generalization 
of the Bhatnagar-Gross-Krook model) and of Anderson- 
Witting @ as special cases if F(uj) — T and F(oj) = ui, 
respectively. 

The stress-energy tensor is 



where 



2T = (?Y) 



corresponds to a perfect fluid and 
IP" = (p>"x) 



(5) 



(6) 



(7) 



encodes dissipative corrections. For a system of massless 
particles we have T M " = e{u»u v + A^/3), where A^ 1 " = 
g^ v + u^u v is the spatial projector. To write down the 
NET j we also introduce a nonequilibriurn tensor 

The 

nonequilibrium tensor is symmetric and transverse, so in 
the rest frame has only spatial components 7^ . Closure 
is obtained from the EPP: 



/ = fa 



1 



T T 2 

r i j , r i j I TTi 
—p P J -f zj + 2T2p2 P P J P V Jijlh 



(8) 



Although the detailed derivation of Eq. ([HJ can be 
found in |23l l2~i| , for the sake of clarity it is convenient 
to briefly sketch the main logical steps that lead to this 
closure. The starting point is the variational equation 
whose solution gives the / that extremizes the entropy 
production 571 [x] with the constraint (imposed by a La- 
grange multiplier) that IP"[x] take on known values. We 
note that Sit and IP" are functionals of x, the deviation 
of / from equilibrium, and that the variation in the vari- 
ational equation is performed with respect to \- The La- 
grange multiplier enforcing the constraint on IP" turns 
out to be proportional to 7^. Solving the variational 
equation (we stay at second order in the relaxation time) 
one finds a relation between x & n d 7^, which, via Eq. 
©, leads directly to Eq. ©. 

Once we have / expressed in terms of our nonequilib- 
rium tensor 7 y , the rest of the calculation is standard: 
we can immediately compute the shear tensor and ob- 
tain the evolution equation for j iJ from the method of 
moments. We leave the latter task for the next subsec- 
tion. 

The first order shear tensor is II 00 = H k0 = 0, and 



IT/ 



The shear viscosity is 



T r T 

~15 



15 



U) 
J2 



(9) 



(10) 



and tv is related by a constant to r r 



G 



U) u 
J2 



The shear tensor at second order is 



3 ' 



Jh 



with 



T 2 
_L / 

45 ^ 



UJ 



(11) 



(12) 



(13) 



It is interesting to remark that if we neglect the second 
term in (JSj) , and set = 0, F = T (corresponding to 
Marie's model) and r r = r)/(e + P), where P = e/3 is the 
pressure, we recover Grad's quadratic ansatz given for a 
Boltzmann gas by 0, Q 



f = fo 1 + 



iwy 

2T 2 (e + P) 



(14) 



This shows that the closure discussed here is a nontrivial 
generalization of the commonly used quadratic ansatz, a 
fact that may have interesting consequences in the study 
of higher order anisotropy coefficients such as the hex- 
adecapole flow V4 (see e.g. [HI). 



B. Evolution equations 

We now turn to the evolution equations of the NET. 

For the hydrodynamic variables (e, u M ) these are sim- 
ply T.^" = 0, i.e. energy-momentum conservation (ex- 
actly as in standard fluid dynamics). The equation of 
motion for 7^ can be found from the transport equation 
([1]) by the method of moments. The second moment of 
the kinetic equation is 



(pW7i/»[7]) = WcoiH) 



(15) 



where we write [7] to emphasize that and I co \ are 
functionals of 7 (recall that they are functionals of \i 
and the solution to the variational equation has provided 
us with a relation between \ an d 7^")- For simplicity, 
in what follows we will neglect the terms cx u^-^j 2 and 
oc ^imlij-,^, which are expected to be very small. We get 

(16) 

where we have defined 



and 



pPp a p~< . . . p K 
pa 



(17) 



(18) 
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for brevity. 

Eq. (|16p determines the evolution of the nonhydrody- 
namic variable 7^. We have already shown in 21] that if 
we expand to second order in veloci ty g radients, the 
above formalism goes over to SOFD [25l - |27| . We empha- 
size that the transport coefficients of the NET are chosen 
to be exactly the same as those of SOFD. This is because 
we want the NET to go over to SOFD when gradients are 
small. In doing so, we are neglecting higher order cor- 
rections to these transport coefficients (see [29|, HH for a 
discussion of this issue) . Since ultimately we wish to ap- 
ply the NET to describe the evolution of strongly coupled 
dense matter, for which kinetic theory does not hold, we 
believe that it is better to regard the transport coeffi- 
cients of the NET as adjustable parameters. In any case, 
once the transport coefficients of SOFD are chosen, the 
NET has no free parameters. 



III. COMPARISON TO KINETIC THEORY 

We now compare the NET to SOFD and to kinetic 
theory. For the sake of simplicity, we focus on 0+1 
boost invariant dynamics, which is a useful toy model 
of a heavy ion collision (see e.g. [H). All quantities are 
independent of space-time rapidity £ = arctanh(z/i) and 
transverse coordinates (x, y) . Proper time is denoted by 

T = Vt 2 -Z 2 . 

It is convenient to introduce a second order transport 
coefficient 



Ai 



45t„-77Z/ 
tvT 3 



(19) 



which arises naturally in the context of SOFD 0, 
The equations of SOFD are (we write II 
a time derivative with an over dot) 



n| and denote 



4 n 

■^-e+~ 
or t 



n 



4t] 



n 4n 



Al r n 2 



(20) 



3 t iT^r] 2 



In the NET, the evolution equation of 7 is (we write 

7 = 7|) 



. = _2 _7 _ 47 _ 2 

St^t T n 3 T T^rj 



(21) 



This equation must be supplemented with the relation 
between II and 7 



II = 2777 + 2A17 2 



(22) 



and with the equation for e given in ([2U)l . We empha- 
size that in the NET one first solves for the nonhydro- 
dynamic variable 7 and then calculates II from it. A 
somewhat similar situation is encountered in the AH of 
[30l - l3"2l | . where the dynamical variable is (apart from the 



hydrodynamic ones) an anisotropy parameter from which 
II is computed a posteriori. 

To solve the Boltzmann equation we use the Lattice 
Boltzmann implementation recently developed in |49j |. 
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Figure 1: (Color online) Evolution of the pressure anisotropy 
Pl/Pt = (Pa - n)/(P + n/2) for 0+1 boost invariant flow 
in the NET, SOFD and kinetic theory (Po is the equilibrium 
pressure), corresponding to the equilibrium initial condition. 
The transport coefficients are set to t„ = 677/ (sT) and Ai = 
r]/(2irT). The upper panel corresponds to rj/s = 0.8, and the 
lower panel corresponds to 77/s = 3.0. 



Figure [T] shows the evolution of the pressure anisotropy 
P L /P T = (P - n)/(P + n/2) computed in the NET, 
SOFD and kinetic theory (Po,l,t are the equilibrium, 
longitudinal and transverse pressures), for the equilib- 
rium initial condition Pl / Pt = 1 ■ The results were 
obtained with the transport coefficients set to tv = 
Qrj/(sT), which is the value obtained from the kinetic 
theory of a Boltzmann gas pl I2H |27| . and \ x = Tj/{2irT). 
We show results for two different values of 77/s, namely 
77/s = 0.8 and 77/s = 3.0. We choose these values be- 
cause they correspond to a typical value for the hadronic 
phase and a very large value used to test our model, re- 
spectively. The initialization time is To = 0.5 fm/c, the 
initial temperature is To = 295 MeV, and the initial en- 
ergy density is eo = 4.9 GeV/fm 3 . Our conclusions do 
not depend on the precise values of these initial values. 

From Figured] we see that the NET can track kinetic 
theory better than SOFD, even for a very large 77/s = 3.0. 
On the other hand, the results obtained with SOFD do 
not agree with those of kinetic theory. For r//s = 3.0, 
the pressure anisotropy obtained from SOFD becomes 
negative, signaling the breakdown of this approach. It is 
worth mentioning that the agreement obtained between 
the NET and kinetic theory for such large values of 77/s 
is similar to that obtained from the developments of [29] 
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and (28|. We note that for low values of 77/s ~ 0.08 (not 
shown in Figure [TJ the agreement between SOFD, the 
NET and kinetic theory is excelent (see Figure [2]) . 
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Figure 2: (Color online) Evolution of the pressure anisotropy 
Pl/Pt = (Po - n)/(P + 11/2) for 0+1 boost invariant flow 
in the NET, SOFD and kinetic theory (Po is the equilib- 
rium pressure) , corresponding to a strong off-equilibrium ini- 
tial condition. The transport coefficients are set to = 
6ri/(sT) and Ai = ri/(2nT). The upper panel corresponds 
to 77/s = 0.08 and 77/s = 0.8, and the lower panel corresponds 
to 77/s = 3.0. 



The equilibrium initial condition used to obtain Figure 
[T] does not correspond to the real situation encountered 
in a heavy ion collision. Matter created in such events 
exhibits very large momentum-space anisotropy at early 
times, so that a highly nonequilibrium initial condition 
for hydrodynamics (and the NET) seems more natural. 
Figure [2] shows the evolution of the pressure anisotropy 
for a strong off-equilibrium initial condition Pl/Pt — 
1.5. This value for Pl/Pt corresponds to |IT | ~ 2P /7, 
which is close to the constraint on the initial Iln , namely 
I IXo j < Po/3, discussed in [l4[ for the reliability of the 
hydrodynamic description of the evolution of the fireball. 

In Figure [5] we also show results for 77/s = 0.08, be- 
cause it is not a priori evident that SOFD and the NET 
will be able to reproduce kinetic theory if a strong off- 
equilibrium initial condition is used. We find that indeed 
this is the case: for such low value of 77/s, the agreement 
between SOFD, the NET and kinetic theory is excelent, 
even for very early times. In contrast, for r)/s = 0.8 the 



SOFD deviates significantly from the result of kinetic 
theory, while the NET still tracks it very well during the 
entire evolution. For the very large value of 77/s = 3.0, 
the pressure anisotropy of SOFD becomes negative for 
t — t > 2 fm/c. On the other hand, the NET is still able 
to track kinetic theory remarkably well at all times. Note, 
however, that slight discrepancies between the NET and 
kinetic theory start to show up. 

As a check, we have used two other set of values for the 
second order transport coefficients, t„ = 3rj/(sT) (strong 
coupling Hi-il) and Ai = t]/{2ttT\, and tv = 677/ (sT) 
and Ai = (standard IS theory [Io|, HU), and found 
that the NET also agrees better with kinetic theory than 
SOFD. 



IV. SUMMARY 

In this paper we have tested the nonequilibrium model 
presented in [23J, |24j against second order fluid dynamics 
with regards to their ability to track the more fundamen- 
tal kinetic theory description of a relativistic system. Our 
model assumes a closure that is nonlinear and does not 
rely on the gradient expansion, thus generalizing Grad's 
quadratic ansatz in a nontrivial way. The evolution equa- 
tions for the macroscopic variables are obtained from the 
kinetic equation by the method of moments. 

By direct comparison to solutions of the Boltzmann 
equation for 0+1 boost invariant flow, we show that the 
effective theory describes the evolution of the pressure 
anisotropy better than standard fluid dynamics. The 
promising comparison to kinetic theory suggests that our 
model could be used to advantage to study the early time 
evolution of matter created in heavy ion collisions, dur- 
ing which very large momentum-space anisotropies oc- 
cur. Moreover, even for a strong off-equilibrium initial 
condition, our model can reproduce the results of kinetic 
theory for very large values of r)/s. These features make 
it an attractive choice to investigate the dynamics of the 
fireball produced in heavy ion collisions in numerical sim- 
ulations including a realistic temperature-dependent 77/s, 
which rises to large values in the hadronic phase before 
freeze out. 
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